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ABSTRACT 

The heating of minor Ions In solar flares by wave- wave-particle 
Interaction with Langmuir waves, or ion acoustic waves, can be described by a 
diffusion equation In velocity-space for the particle distribution function. 
The dependence of the heating on the Ion charge and mass, and on the 
composition of the plasma, is examined In detail. It Is found that th^ 
heating mechanisms proposed by Ibragimov and Kocharov (1977) and Kocharov and 
Kocharov (1981) cannot account for the enhanced abundances of heavy elements 
In the solar cosmic rays. 

Subject Headings: cosmic rays; general -- particle acceleration — 

turbulence 
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I. INTRODUCTION 

The heating of plasma particles by turbulence Is a consequence of the 
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dissipation of wave energy. Often, wave damping results from the resonant 
Interaction of a small population of particles with the waves, such as In the 
case of Landau damping of electrostatic waves. This type of damping affects 
those particles with velocities matching the phase velocity of the wave, I.e., 
a wave with frequency w and wavenumber k will Interact with a particle of 
velocity V for which w - k«v » 0. Because of their small thermal velocities, 
positively charged Ions contribute little to this If.teractlon. Another kind 
of wave-particle Interaction takes two waves, and can directly Involve thermal 
Ions. The resonance condition for the particles and waves of frequencies w 
and w* and wavevectors k and k’ respectively, Is cu - k*v * oi' - k'*v. 

The loss of wave energy v, this process Is referred to as "nonlinear" Landau 
damping (Davidson, 1972). 

The theory of nonlinear Landau damping has recently been applied to the 
cosmic ray problem. In particular to the question of Ion abundances In solar 
flare events. The generalization of the theory to a plasma containing man' 

1on species Is made In this paper. 

It has been known for many years that cosmic rays from the sun have a 
composition which Is different from typical solar abundances (McGuire, 
Rosenvinge and McDonald, 197d; Mason, Fisk, Hovestadt and Gloeckler, 1980). 

For example, energetic particles (In the MeV range) will ocasslonally show an 
extraordinary enrichment of ^He. The ratio ^He/'+He has been observed as large 
has been observed as large as unity, which Is an enhancement of ^He of 10*+ 
over Its normal solar abundance. Increased abundances of the heavy elements 
seem to be characteristic of flares, and Is well documented for nuclei up to 
atomic number twenty-six. However, the relative abundances of these elements 



are different from flare to flare, while the greatest enhancements are 
observed In flares that are not large. 
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Attempts to explain different abundances based on nuclear processes 
occuring In flares run Into difficulty because of the lack of deuterium or 
tritium by-products. Several theories were made for a selective acceleration 
mechanisms (Ibragimov and Kocharov, 1977; Fisk, 1978; Kocharov and Kocharov, 
1978, 1981; Hayakawa, 1983; Varvoglis and Papadopoulos, 1983). 

There Is some Interest among the cosmic ray community In the acceleration 
mechanism proposed by Ibragimov and Kocharov (1977) and Kocharov and Kocharov 
(1978 and 1981), which are based on stochastic acceleration by Langmuir and 
Ion acoustic waves. However, on close examination, one finds that the 
formulas which these authors use to describe the turbulent heating do not 
contain the proper dependence on Ion charge.. For this and other reasons to be 
discussed In this paper, this heating mechanism cannot lead to enrichment of 
^He and other heavy elements. 

The turbulent heating mechanism can be viewed as a scattering process. A 
charged particle, set Into oscillatory motion by a plasma wave, will radiate 
Into other plasma modes, and absorb the residual wave energy and momentum. In 
obvious analogy with the conventional theory of the scattering of 
electromagnetic radiation on charges, this Is known as plasma Compton 
scattering. Furthermore, because the charged particle displaces surrounding 
charges to produce a screening cloud, there Is also the possibility of 
scattering on these charges. This Is known as "nonlinear scattering". These 
mechanisms are reviewed below. 


II. THEORY OF SCATTERING 

There are two methods to study the scattering problem which lead to 
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Identical results. One Is the tot particle method which, as Its name would 
suggest, solves for the motion and radiation due to an Individual '’sample” 
particle. Wave-particle Interactions have been studies extensively with this 
theory based upon a semi -classical formalism for the scattering processes 
(see, for example, the books by Tsytovich, 1970 and Melrose, 1980). A second 
approach uses a perturbation theory which Is based on the Vlasov cumulant 
hierarchy (Davidson, 1972). 

The test particle method can be used to calculate the probability of 
scatterlno a wave with wavevector k Into a wave with wavevector k* by a 
particle of charge Q and mass M, moving with velocity v. The derivation can 
be found In Tsytovich (1970) or Melrose (1980), but basically It Involviis 
calculating the work done on the field E (k‘) by the currents which arise from 
the forced oscillation of the charge In the field E (k). It will suffice to 
remark that there are two contributions to the current. There Is the current 
due to the particle Itself. To this It Is necessary to add the current due to 
the polarization response of the plasma to the test charge. The current due 
to the oscillation of the test charge leads to Compton scattering (see Figure 
1), while the current associated with the screening charges around the test 
charge produces the nonlinear scattering (Figure 2). Nonlinear scattering can 
Interfere destructively with the Compton scattering since the currents add 
linearly. 

The calculation of the scattering probability forms the basis of the 
theory of wave scattering and particle heating. For this theory, it is useful 
to define the wave occupation number 
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E(k) 1s the Fourier transform of the field E(r,t); 


E(k) • / e E(r,t) d^r. (2) 

In the semi -classical language, Hw is the energy of a quantum of wavenumber 
k, i\ being Planck's constant divided by 2tt. e{o),k) Is the dielectric function 
of the plasma. (The normal modes of the plasma occur for those frequencies 
o)(k) for which the dielectric function gives zero.) The volume V of the 
plasma Is taken In the Infinite limit. As defined here, Nj^ Is 
dimensionless. In terms of Nj^ , the total energy In the waves (electric and 
kinetic) Is given by 

W » V / -ftu(lt) (3) 


The probability for scattering a wavevector k Into a wavevector k‘ by a 
particle with momentum p Is the following (Melrose, 1980); 


oi(p. k, kM * 4(2,t)3 Q2 [ - - fe" ' 4(a.-(.'-(k-k‘).¥). 


5aT 




(4) 


Melrose^ gives the scattering matrix element A for the case of a hydrogen 


^See Melrose (1980), eq. (10.114). The factor qm/emg in this equation should 
be em/qm^: this can be verified by adding of eq. (10.112) with 

of eq. (3.96), which correspond to the currents associated with the 
polarization cloud and the test charge, respectively. 


plasma In which electrons dominate the nonlinear scattering. In a solar 
environment, the plasma consists mainly of protons, electrons, and helium (the 
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ratio of the number of helium nuclei to protons Is the parameter n^^/n^ * a). 
To study the heating of minor constituent species, It will be necessary to 
include the nonlinear scp.ttering by the positl^'ely charged ions In the 
screening cloud. It Is useful to write the matrix element as the sum of two 
parts; 

A * Ag + Agp. (5) 

The first term comes from Compton scattering and Is given by (3.99) of Melrose 
(1980); 
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The terms of order kv/oj and k'v/o)' are small, and are called "Doppler" 
corrections to Compton scattering. The term Agp derives from the nonlinear 
scattering. To evaluate this term. It will be assumed j-n1y Ion species 

with thermal speeds less than the phase speeds of the waves contribute 
significantly to the scattering. It Is also taken for granted that the speed 
of the test particle Is less than the wave speed. The nonlinear scattering 
matrix element Is then found to be 



k*k' 1 


I j—Xy (k-k', oi-o)'). 
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V denotes the sum over the Ion species for which the thermal velocity Vj^ is 
less than oj/k and w'/k' Nots that the subscripts on charge q, mass m, number 
density n, temperature T, and other quantities denote the plasma species; e 1s 
used for the elementary charge, e = Iq^l. The thermal speed of species y Is 
defined (KT^/m^)^'^^ , where k is Boltzmann's constant. The plasma 
susceptibility of species y Is given by 


x^(k.w) 
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where f Is the normalized velocity distribution function. The dielectric 
function can be written as a sum over all the Ion susceptibilities; 


elk.o)) *1 + 1 X„ (kfW). (9) 

u ^ 

The plasma physics can be simplified when the group velocity, 
o)-o)7|lc-lc' I , associated with the scattering Is less than the thermal velocity 
of tha plasma particles. This Is true for most scattering Ions because the 
resonance condition embodied In the delta-function In equation (4) requires 
that the group velocity be less than the velocity v of the test charge. 
Therefore, as long as v Is not too large. It follows that 

a^rn 0 ^ 

x„(ic-i(‘. ^-0.') [ v, - fcV — rf-" 

for the constituent plasma particles. 

The rate of change of the particle distribution function due to the 
scattering process can be written In terms of the scattering probability In 
the form of a diffusion equation (Melrose, 1980): 
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Ibragimov and Kocharov (1977) Introduced the Idea that this process could 
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be a mechanism by which unusual particle abundances could be produced 1n solar 
flares. Their theory presupposes an acceleration process which Is very 
efficient for particles above a certain threshold velocity (such as the sound 
speed or some other hydrodynamic speed). The role of the nonlinear scattering 
Is to heat the Ions, thereby enhancing the number of thermal particles with 
speeds larger than the threshold speed. Ibragimov and Kocharov considered 
heating caused by Langmuir waves. Later work by Kocharov and Kocharov (1978, 
1981) applied the theory to Ion acoustic turbulence. They also Include linear 
scattering as a secondary acceleration process. 

The basis of these theories Is the diffusion In velocity of the 
distribution F of plasma test particles. Consider the diffusion equation (11) 
In spherical polar coordinates, assuming Isotropic distribution F and M|^: 



For constant D, there Is a simple solution for an Initial Maxwellian 
distribution given by 

F = exp(-vVvg) ; (14) 

v2 . ^ + 4Dt . (15) 

The quantities T, M, and D are associated with a particular species of test 
particle. It Is then an easy matter to calculate the number of particles 
Injected Into velocities greater than the critical velocity, v^: 

00 

n* = n / F(v) 4*v2dv . 


(16) 


Naturally, particles with contparatively large diffusion coefficients will 
experience a qreater enrichment In n*. The enrichment for different elements 
depends upon D, which Is a function of char^^e and mass. 


III. SCATTERING ON LANGMUIR WAVES 

The first case which will be considered Is the scattering of Langmuir 
waves on charged particles. Langmuir waves are longitudinal plasma waves 
having the following dispersion relation; 


0)2 ■ + 31(2 , (X7) 

The electron plasma frequency Is o» » (Atrn.e^/m.)^^^. This mode occurs only for 

0 0 0 

wavenumbers k much less than the Inverse Debye length, ^ (Awnge^/KT^)^'^^. 
Langmuir waves have relatively fast phase velocities; 
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Generally, the phase speed 1s larger than any thermal speed In the plasma. 

It Is well known that for scattering on electrons, Compton scattering and 
nonlinear scattering of Langmuir waves Interfere destructively (to lowest 
order In kv/o)^). This can be deduced from equation (5). 

In the case of a heavy Ion with charge Q and mass M, the nonlinear 
scattering Is much larger than the Compton term. This Is because the 
nonlinear scattering takes place on the lower mass electrons which screen the 
Ion, rather than on the ion Itself. The scattering probability from equation 
(4), together with (10) for the particle susceptibilities, gives for Langmuir 


waves 
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Use Is made of the fact that 3c(k,u)/3w » 2/w^ for Langmuir waves. The result 
also follows from the derivation leading up to (7.71) of Tsytovlch (1970). 

The diffusion coefficient for scattering on Langmuir waves can be 
Inferred from equation (12) and the scattering probability given above. Since 
nothing Is known about the Langmuir distribution, It Is assumed that waves 
fill k-space uniformly for wavenumbers less than some minimum wavenumber 
(which must be less than As Iona as the Ion velocity Is less than 

(kXg)vTe» the diffusion coefficient Is 


0^ » .16 tt 2 Z2e2 (W^)2/(A2mVn^KT^(l + Te/T„)2). (20) 

P 6 6 6 P 

Is the energy density In Langmuir waves as defined In eq. (7). Z Is the 
charge ratio Q/e, and A the mass ratio M/mp of the scattering Ion to the 
proton. The factor .16 derives from the Integration over the distribution of 
waves, which of course depends upon the assumed spectrum. The diffusion 
coefficient defined here Is not proportional to Z^/A^ as In Ibragimov and 
Kocharov (1977), but rather Z^/A^. The formula of Ibragimov and Kocharov Is 
presumed to be In error. The charge to mass factor derived here does not 
favor the preferential heating of the heavier Ions. 

Another drawback Is related to particle collisions. If the effect of the 
turbulent scattering Is to heat Ions differently, the effect of collisions 
between Ions will be to keep them at equal temperatures. In order to compare 
these competing processes, one can compute the diffusion time scale and the 
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collision tfnie scale, The ambient plasma Is assumed to correspond to the 
upper chromosphere In which electrons have been rapidly heated (Ibragimov 
andKocharov, 1977, and Kocharov and Kocharov, 1981). The density and 
temperature are given by n^ « 10® cm"®, Tp » SxlO** K, and T^ » 100 Tp. A 
speculative value (Hoyng, 1977) for the turbulent energy In Langmuir waves 
during solar flares can be W^/ngicTg * 10 “ 2 , For these parameters, the 
diffusion rate Is 




sec’^ 


( 21 ) 


Therefore, the Ion distribution will heat up within about 16 seconds. The 
time for equipartitlon (Spitzer, 1962) between helium and protons Is teq » 
6xl0“3s fo this plasma* As long as the equipartitlon time Is so much faster 
than the heating, It Is hard to see how any appreciable temperature difference 
between Ions can occur from the nonlinear scattering. The collisions would be 
less Important If the Ions were at a temperature 10® K or larger; but In this 
case the electron temperature would also have to be larger In order to keep 
the Injection velocity (the sound speed) for the secondary acceleration 
mechanism well above the Ion thermal velocity. The low Ion temperature Is 
also preferred In this theory to account for the low Ionization of particular 
Ions which are not greatly enhanced In these events (Kocharov and Kocharov, 
1981). 


IV. SCATTERING ON ION-ACOUSTIC WAVES 
Ion acoustic waves can also scatter the Ions, This wave mode has the 


dispersion relation 
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u) * kCg 


( 22 ) 


where the sound speed Is approximately Cj « (Tg/mp)V^» Sound waves are 
heavily damped unless Tg » Tp. It follows that when this wave exists 1t has 
a phase velocity which Is much larger than the thermal velocity of protons 
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However, the therma'J velocity of electrons Is greater then the phase speed. 

For this reason the nonlinear scattering Is predominantly due to positive Ions 
In the screening cloud. 

One again, the group velocity of the scattering waves Is less than the 

thermal velocities. Therefore, the value for the susceptibilities at the beat 

frequency /ji-u Us given in the limit of small frequency by equation (10). 

For th>ji case when the scattering Ion is a proton In a pure hydrogen 

plasma, It Is seen from equation (5) that the screening by the Ions exactly 

cancels the Compton scattering to zero order In kv/w. Note, however, that the 

cancelation does not occur for heavier Ions which have charge to mass ratio 

different from protons. Moreover, the cancelation will not occur If there Is 

significant screening by other Ions In addition to protons. 

From the Ion acoustic wave dispersion relation and equation (6), one can 

see that the correction to Compton scattering due to the "Doppler effect" Is 

approximately v/Cg, which Is generally small for heavy thermal Ions, The 

corrections to the screening term due to the presence of other positively 

charged Ions are of order n /n_, the ratio of the number density of ton u to 

P 

protons. Corrections due to the screening by electrons are of order Tp/Tg. 

The scattering probability for Ion acoustic waves 1s the following: 
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In deriving this fornula, terms of order Tp/T^ have been negelected. The 
relation 3c(k,o))/3u) » Sirnpe^/rfipW^ for Ion acoustic waves has also been used. 

The diffusion coefficient for Ion acoustic waves follows from the 
scattering probability. The diffusion coefficient will be split Into two 
parts D® » Dj + Djj. Assuming an Isotropic wave spectrum, 


Dy » .90 tt 2 C / k® Nj dk; 

Vgr K 


(25) 
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The group velocity of Ion acoustic waves is The factor C Is defined 


r .4^2 Z‘^e‘* 1 

^ ]^(7irV Pn?rre2nr)2 • '2/) 

P P 

The diffusion coefficient Dj was written down by Kocharov and Kocharov 
(1978), except that their factor C Is proportional to zVa^, Instead of the 
factor Z^/A^ given above. There Is no physical reason for this difference. 
Dj Is due to the Doppler correction to Compton scattering. 

The term Djj arises from the Incomplete destructive Interference of the 
scattering from the screening Ions with the Compton scattering. When the 
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Helium component In the plasma Is disregarded, Djj Is zero for Ions Z/A»l. 

For this reason, it may have been overlooked In the Kocharov and Kocharov 
theory. In a mixed plasma, Dji Is generally larger than Dj for Ions with 
velocities less than thermal. 

The Integral over occupation number Is needed to evaluate the diffusion 
coefficient. For this, the spectrum of Ion acoustic waves Is presumed to be 
peaked around some wavenumber k^. For simplicity, Is taken to be constant 
between k^ + and k^ - where Ak Is the width of the spectrum. The 
diffusion coefficient Is 
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where 
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It Is assumed that the spectrum Is not too narrow: Ak/k^ should be larger 

than v/Cj. 

To facilitate comparison with the heating by Langmuir waves (see equation 
21), the diffusion coefficient Is evaluated for the same parameters. In 
addition W^/n^xT^ = 10 "^, k^x^ » Ak/k^ =1/^ , and o * 20 %: 

nS 

8.1 X 103 5 sec 1. . (30) 

D Is plotted for three Ion species In figure 3. The diffusion is much more 
efficient than In the Langmuir scattering case, and the diffusion time is even 
less than the time for Ion-Ion collisions for the level of turbulence. 

However, this acceleration mechanism still has fundamental problems. It can 


be seen from figure 3 that the bulk of ^He Is actually scattered less than ^He 
or hydrogen. This Is due to screening effects by the ambient helium, which 
makes up 20% of the plasma. It Is also apparent that the scattering of the 
heavy elements Is not enhanced by the charge and mass dependence of the 
scattering probability. 


V. SUMMARY 

The enchanced abundance of the heavier Ions In solar flares had been 
attributed by Ibragmov and Kocharov (1977) and Kocharov and Kocharov (1981) to 
particle scattering on plasma turbulence with an efficiency proportional to 
Z^/A^. However, as discussed In this paper, the diffusion coefficient 
associated with Compton scattering Is proportional to 

® 

while the diffusion coefficient for nonlinear scattering Is proportional to 

Z2 7*2 

where V and A' are the charge and mass numbers Of the plasma Ions In the 
charge cloud around the scattering particle. Therefore, the dependence 

has no basis In this theory, and Its use probably derives from an Incorrect 
application of Tsytovich's theory (1970) to multiply charged, heavy Ions. 

This point is crucial for the application to solar cosmic ray theory, since 
without the Z‘*/a 2 dependence, the nonlinear scattering on ion acoustic or 
Langmuir waves cannot account for the enhancement of the heavier Ions. 

In order to produce any significant temperature difff between Ion 



species, the rate of heating by the turbulence must be treater than the rate 
at which collisions between Ions tend to equilibrate the temperature. For the 
parameters assumed to pertain to solar flares, together with estimates for the 
turbulent energy and spectrum, this does not seem to be true for the nonlinear 
scattering on Langmuir waves. 

On the other hand, the scattering by Ion acoustic waves seems to be much 
more efficient than scattering by Langmuir waves. In a mixed Ion plasma, 
however, the constructive Interference between the Compton and nonlinear 
scattering on Ion acoustic waves Is sensitive to the ratio of Z/A. For the 
solar plasma, this effect turns out to be detrimental to the heating of ^He. 

In conclusion, the acceleration mechanism proposed by Ibragimov 
andKocharov (1977) and Kocharov and Kocharov (1978, 1981) does not have the 
sensitivity to ion charge needed to account for the observed solar cosmic ray 
abundances* 
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Figure Captions 


Figure 1. Conpton scattering on the Ion Q of a plasna wave with 
wavenumber k Into a wave with wavenumber k’. 

Figure 2. Nonlinear scattering of a plasma wave with wavenumber k Into a 
wave with wavenumber k* through Intereactlon with the charge cloud of Ion 
Q. 


Figure 3. The factor D, which appears In the diffusion coefficient (equation 
29) associated with scattering on Ion acoustic waves* as a function of 
velocity for three Ion species. The helium ratio Is a * 0.2. Velocity Is 
In units of the proton thermal velocity. The sound speed Is Cj = 10 vjp. 
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